Magnetism and Electronic Correlations in Quasi-One-Dimensional Compounds 
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In this contribution on the celebration of the 80th birthday anniversary of Prof. Ricardo Fer- 
reira, we present a brief survey on the magnetism of quasi-one-dimensional compounds. This has 
been a research area of intense activity particularly since the first experimental announcements of 
magnetism in organic and organometallic polymers in the mid 80's. We review experimental and 
theoretical achievements on the field, featuring chain systems of correlated electrons in a special 
AB2 unit cell structure present in inorganic and organic compounds. 
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I. INTRODUCTION 

The magnetism of organie^'^'^ and organometallio^ 
polymers has been a challenging topical field since its 
first experimental announcements. This rapidly growing 
and interdisciplinary research area also includes inorganic 
compounds- , with ferro- and ferrimagnetic long-range or- 
der (see Sec. II below). In this work we briefly review 
some attempts to describe the ground state and the low- 
temperature thermodynamics of these compounds. In 
particular, we report on analytical and numerical results 
on polymeric chains of correlated electrons in special unit 
cell topologies shown in Fig. [TJ 



II. QUASI-ONE-DIMENSIONAL MAGNETIC 
COMPOUNDS: A BRIEF REVIEW 



Despite many years of experimental and theoretical ef- 
forts, the complete understanding and precise charac- 
terization of magnetism and electronic correlations in 
quasi-one-dimensional (quasi-ld) compounds still offer 
great scientific challenges and technical difficulties^. Re- 
garding, for instance, organic magnetic polymers, it is 
known^ that their magnetic properties are ascribed to 
the correlated p-electrons of light elements, such as C, 
O, and N, in contrast to the magnetism found in tran- 
sition and rare-earth metals due to partially filled d or 
f orbitals. For this reason it took many years of efforts 
on the synthesis and characterization of a great variety 
of compounds before the announcement of bulk ferro- 
magnetism (FM) in an organic polymer—. In Fig. \V[h) 
we sketch this material made of polyacetilene-based radi- 
cals, R* , containing unpaired residual electrons, i.e., poly- 
BIPO or poly[l,4-bis(2,2,6,6-tetramethyl-4-piperidyl-l- 
oxyl)-butadiyne] . However, this compound presented 
several problems due to its insolubility and poor repro- 



ducibility both in the preparation and in the magnetic 
results^. Later, Nishide and collaborators^ have suc- 
cessfully synthesized polyphenyacetylenes with various 
types of radical groups. These polymers exhibit sim- 
ilar band structure schemes^ comprising filled bonding 
molecular-orbital bands, empty antibonding bands, and 
narrow half-filled nonbonding bands, usually just one at 
the center of the band. A net magnetic moment may 
appear either because the number of itinerant antiferro- 
magnetically (AF)-correlated w electrons per unit cell is 
odd and/or due to the presence of localized electrons-'^-. 

A seminal work has also been performed by Takahashi 
and collaborators in order to extensively characterize the 
long-range macroscopic FM behavior found in the or- 
ganic compound p-nitrophenyl nitroxyl nitroxide radi- 
cal (p-NPNN, in the 7 and (3 phases)^. Actually, the 
excellent fitting of the low-temperature (T) experimen- 
tal data is consistent with predictions from the thermo- 
dynamic Jiethe- ansatz solution of the Id-quantum FM 
Heisenberg model^: susceptibility x ~ T~^ and specific 
heat C - Ti/2, as T ^ 0. 

Other organic magnets have been synthesized, such as 
polyradicals derived from poly(l,3-phenylenemethylene) 
and polyphenylenevinylene-based radicals'^. In these 
cases the polymer structure is made of benzene rings 
linked by divalent carbon atoms or including pendant 
radicals with oxygen atoms carrying an uncompensated 
electror*^. Another family of organic magnetic polymers 
is that of the doped poly(TO-aniline) compounds, in which 
the carbon atoms responsible for the links between the 
benzene rings are substituted by ionized nitrogen with 
a H-bond or a radical plus a charge acceptor—. On the 
other hand, doped polypyrrole compounds also exhibitii^ 
interesting magnetic properties and Drude metallic re- 
sponse as well. 

A distinct class of magnetic polymers combines 
metal ions with organic complexes, displaying a rich 



variety of magnetic behaviors, such as ferro- and 
ferrimagnetism, AF and canted AF and spin glass 
phaseiii^. In fact, the first experimental observa- 
tion of a magnet with spin residing in a p-orbital 
was performed in the compound [Fe(C5Me5)2]^[TCNE]^ 
(TCNE ~ tetracyanoethylene)^. Some homometallic 
fcrrimagnets with chain structure^'^'^'^ involve the com- 
poundsi5-i6 M2(EDTA)(H20)4-2H20 {M = Ni, Co; 
EDTA = ethylenediamminetetraacetate — C10N2O8) 
and M(i?-py)2(N3)2 (M =Cu, Mn; R-py = pyri- 
dinic ligand = C5H4N-i? with R = CI, CH3, 
etc.)^>^i^'^. Regarding bimetallic chain materials, 
the compounds! MnCu(pbaOH)(H20)3 [pbaOH = 2- 
hydroxy-1, 3-propylenebis(oxamato) — C7H6N2O7] has 
been one of the first synthesized which retains long- 
range FM or ferrimagnetic order on the scale of the 
crystal lattice, as in the case of isomorphous realiza- 
tions^^iSSii^. Heterometallic chain structures have also 
been object of systematic study^''"-^-. More recently, the 
metal-radical hybrid strategy, combined with fabrication 
of novel polyradicals^, has led to the synthesis of a vari- 
ety of heterospin chain compound a^^'^^ . Several of these 
compounds display Id ferrimagnetic behavior-^Si^i mod- 
eled by alternating spin chains^, such as, for instance, 
those with the structure shown in Fig. [TJc) . 

Of recently growing interest we mention the quasi-ld 
chains with AB2 and ABB' unit cell structure [hence- 
forth referred to as AB2 chains; see Figs. [T{a) and (b)]. 
Such structures are found both in inorganic and organic 
ferrimagnetic compounds. Regarding inorganic materi- 
als, we cite the homometallic compounds with a line of 
trimer clusters characteristic of phosphates of formula 
A3Cu3(P04)4, where A = Ca,^'^^'^, Sr,^>^^>S, 
and pb 35,36,38 ^ r^Yie trimers have three Cu"'"^ param- 
agnetic ions of spin S = 1/2 AF coupled. Although 
the superexchange interaction is much weaker than the 
intratrimer coupling, it proves sufficient to turn them 
into bulk fcrrimagnets. Furthermore, compounds of for- 
mula Ca3_2;Srj;Cu3(P04)4, < a; < S'^^, hybrid analo- 
gous to the mentioned phosphates, have also been syn- 
thesized in an attempt to tune the AF bridges be- 
tween Cu sites and possibly explore how paramagnetic 
spins grow into bulk ferrimagnets. It is also inter- 
esting to mention the frustrated AB2 inorganic com- 
pound'^- Cu3 (003)2 (011)2, which displays low-T short- 
range magnetic order and has its physical properties 
well described through the distorted diamond chain 
model4^. At last, we observe that the AB2 struc- 
ture is also presenli^ in the organic ferrimagnetic com- 
pound 2-[3',5'-bis(Af-iert-butylaminoxyl)phenyl]-4,4,5,5- 
tetramethyl-4,5-dihydro-li/-imidazol-l-oxyl 3-oxide, or 
PNNBNO, consisting of three S = 1/2 paramagnetic rad- 
icals in its unit cell. 
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FIG. 1: Ferrimagnetic ground-state configurations of (a) bi- 
partite lozenge AB2 chains, (b) substituted polyacetilene, 
with lateral radicals R* as B sites containing unpaired resid- 
ual electrons, and (c) alternate spin chains. 



III. MAGNETIC CHAINS WITH AB-2 UNIT 

CELL TOPOLOGY: ANALYTICAL AND 

NUMERICAL STUDIES 



In this section we model and discuss analytical and 
numerical results on the magnetic chains with AB2 and 
ABB' unit cell topologies displayed in Figs.[TJa) and (b). 
Eventually, alternate spin chains shown in Fig. [BJc) are 
also considered. 

A rigorous theorem by Lieb^ predicts that bipar- 
tite AB2 chains modeled through a Hubbard Hamil- 
tonian [see Eq. ([T|) below], with one electron per site 
on average (half-filled limit) and repulsive Coulombian 
interaction, present average ground-state spin per unit 
cell h/2 and quantum ferrimagnetic long-range order at 
T = 42,43,44,45 _ fpj^g magnetic excitations on this state 
have been studied in detail both in the weak- and strong- 
coupling limits^, and in the light of the quantum AB2 
Heisenberg model^i^i^. Further studies have consid- 
ered the anisotropic^ and isotropic^*' critical behavior of 
the quantum AB2 Heisenberg model, including its spher- 
ical version'^ ^, and the statistical mechanics of the classi- 
cal AB2 Heisenberg model^. 

Away from half-filling, doped AB2 Hubbard chains 
were previously studied through Hartree-Fock, exact di- 
agonalization and quantum Monte Carlo techniques both 
in the weak- and strong-coupling limitg^'^'-^, including 
also the t — J modelS using the density-matrix renor- 
malization group and recurrent variational Ansdtzes, and 
the infinite Coulombian repulsion limilj^ using exact di- 
agonalization. In particular, these chains represent an 
alternative route to reaching 2d quantum physics from 
Id systems^'^. 



A. Analytical Results 

We start considering the AB2 chain modeUed through 
the one-band Hubbard model, which is the simplest lat- 
tice model for strongly correlated materials: 
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where c]^ ^ (cia^a) is the creation (annihilation) operator 
for electrons with spin (j{=1, |) at site a = A, Bi or B2 
of the unit cell i, i is the hopping parameter, U is the 
intrasite Coulomb repulsion and, in the first summation, 
ia and jP are nearest neighbor sites. We define N as the 
total number of sites and Nc{~ N/3) as the number of 
unit cells. We remark that in the limit U = 00^, the 
Hubbard Hamiltonian reduces to a hopping term with 
double site occupancy excluded. 

In the tight-binding description {U — 0), this model 
presents three bands: one flat with Nc localized or- 
bitals with energy e = and two dispersive with e± = 
±2V2cos{q/2), where q = 2ttI/Nc and I = 0, ...,Nc - 1. 
At half-filling (A^e = N, where A^e is the number of elec- 
trons) and U = the ground state (GS) total spin quan- 
tum number Sg is degenerate, with 5*^ ranging from the 
minimum value (0 or 1/2) to Sg = \Nb ~ Na\/2 = Nc/2, 
where Na {Nb) is the number of sites in the A (B) sub- 
lattice. As proved by Lieb*^^, in the general case of any 
bipartite lattice with Nb ^ ^Va, the Coulomb repul- 
sion lifts this huge degeneracy and selects the state with 
Sg = \Nb — Na\/2 to be the unique GS of the system, 
apart from the trivial {2Sg + l)-fold rotational degener- 
acy. 

In the strong coupling limit ([/ >> t) and at half- filling 
the AB2 Hubbard Hamiltonian, Eq. ([T]), is mapped^i^ 
onto the quantum 5 = 1/2 Heisenberg model with 0(n), 
n = 3, rotational symmetry: 
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where the localized spins Sm interact antiferromagnet- 
ically through Jff = J = U"^ /U > 0. In fact, Lieb 
and Mattis have shown'^^ that the Heisenberg model in a 
bipartite lattice has also Sg — \Nb — Na\/2, which indi- 
cates that Sg — Nc/2 for the AB2 Heisenberg model, 
as in the Hubbard case. Eventually, in the presence 
of an uniform magnetic field H along the z direction 
a Zeeman energy term, —giJi'sH/h^^^Sfc,, is added to 
Eq. (12) , where g is the gyromagnetic factor and jib is the 
Bohr magneton (in what follows we take units in which 
g^B = !)• In 7f = the ground-state of the system 
exhibits^i^ unsaturated FM or ferrimagnetic configura- 
tions as indicated in Figs.[T](a) and (b), with average spin 
per unit cell {S^^n) = fi/2 (Lieb's thcorcm^^). In addi- 
tion, we can also model the alternate spin chains shown in 
Fig-IHc) by considering in Eq. ^ SiQ = Si and S.;^ — s^, 
with S > s. Such model systems have been used to 



describe^i^ a number of organometallic compounds in 
which, e.g., 5 = 5/2 or 5 = 2 and s = 1/2. 

The Euclidean action of the partition function 
in a coherent-state n-ficld representation's., Z = 
J'Dnexp{-SE/fi), is given by Se = S^xc + Swz + Sz, 
with contributions from exchange and Zeeman interac- 
tions, and a topological Wess-Zumino term, which is a 
Berry's phase-like term associated with the time evolu- 
tion of the spin due to quantum fluctuationSi^. The low- 
lying properties of the quantum ferrimagnetic chains in 
Fig. [1] are dominated by infrared fluctuations around the 
Neel configuration. In order to obtain their effective low- 
lying action, we take staggered dimensionless unit coher- 
ent magnetization fields and split the topological term 
in a FM and an AF contribution, Swz = S^^z + ^w^L 
following the spin structure of the polymer. Then, tak- 
ing the continuum limit and integrating out the rapidly 
fluctuating field modes, we obtain the low-energy ef- 
fective action, Seff = L dx/ {2a) L drC, where L is 
the length of the chain of lattice parameter 2a, and 
(3 = {kBT)~^ = it/h expresses the result of a Wick ro- 
tation to imaginary times it. The Lagrangian density is 

C = CaNL + -CS,^! + Cz + CirreU whcrC 
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corresponds to the quantum nonlinear (NL) a model, 
with unit magnetization fields m^ = 1, 
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represents the contribution from the FM Wess-Zumino 
term, and Cirrei are irrelevant terms in the renormal- 
ization group (RG) context. In Eqs. ^ and ^, ax — 
2JS'^H^a^, ar = 1/(8J), and S = S for the AB2 chains, 
whereas ax = JS^tfa^, ar = s/{AJS), and S = S — sior 
the alternate spin chains of Fig.[T] The FM Wess-Zumino 
term is responsible for the ferrimagnetic ground states. 
Indeed, C^J^^ = in either cases of AB or S — s chains. 
Those represent usual quantum AF Heisenberg chains, 
which are known to follow Haldane's conjectured, i.e., 
half-integer spin chains are critical with no long-range 
order, and integer spin chains are disordered. Moreover, 
the addition of the relevant Wess-Zumino term to the 
quantum NL a model, Eq. ([3]), changes its properties 
dramatically since the critical dynamical exponent as- 
sumes z = 2 (nonrelativistic feature), in contrast with 
the value z = 1 found in the relativistic quantum NL a 
model, associated with the 2d-quantum AF Heisenberg 
Hamiltonian. In fact, Eq. ^ corresponds to the field- 
theoretical version of the topological constraints imposed 
by the polymer structure, as identified by semi-empirical 
mcthodsii^°-6^ 

We now perform a momentum-shell low-T RG study of 
the system (see^ and^ to similar treatments to the clas- 
sical and quantum z = 1 NL a models). First, we decom- 
pose the magnetization fields into transversal and longi- 
tudinal components, integrate over the latter one, expand 



the resulting action, Seff = 5*^^^ + S^'^^ + ..., and Fourier 
transform the terms to the momentum-k and Matsubara 
frequency-Wn space, with ujn — 27rn/u, n = 0,±1,±2..., 



u = CP, and C = ir'^JSh? or C = n'^ J S sh^ / {S - s) respec- 
tively for the AB2 or alternate spin chains. The quadratic 
term in the diagonal field space {(/)*, 0} reads^ 
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where hg^ = H(/h, the density of degrees of freedom p 
comes from the integration over the longitudinal compo- 
nents, and the meaning of A^ is discussed below. The 



bare coupling is defined as go = t^/S in rf = 1. The 
quartic contribution is given by^ 
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In the sequence, we require that the fluctuation modes 
in the two-point vertex function scales homogeneously 
through a RG scaling transformation. A; — > 5fc, a;„ — > 
&^a;„, with b = e^. We also take the fixed point AjJ = 
for d = 1, as a consequence of the irrelevance of the uj'^- 
dependent terms in the RG context. The one-loop equa- 
tions for the renormalized coupling g, and dimensionless 
temperature t = g/u and magnetic field h = hg m d = 1 
and z = 2 read^: 



dg_ 
di 
dt 

di 
dh 



-{d+z-2)g+ ^g^ coth[g(l + h)/2t] , 
-(d-2)f+^5tcoth[5(l + /i)/2i] , (7) 



2h 



du 



-zu , u 



We thus obtain the semiclassical fixed point: g* = t* = 
h* — 0, since g* = implies in 5 — *■ 00, and the quan- 
tum critical fixed point^: g* = gc ~ 2tt, t* ^ h* = 0. 
The former describes a Id-classical Heisenberg ferromag- 
net with quantum corrections, whereas the latter is iden- 
tified with a classical Heisenberg model in d -I- z = 3 
dimensions. The analysis of stability shows that both 



fixed points are unstable under thermal fluctuations, but 
only the semiclassical fixed point is stable under infrared 
quantum fluctuations, as shown in Fig. [21 

By studying the correlation length ^ and magnetic 
susceptibility x "^^ identify^S three distinct quantum 
regimes. As T ^ and g > g^, the quantum z = 2 
NL a model is in a quantum disordered phase, whereas 
for g < gc its ground state has long-range order, with 
both quantum and thermal fluctuations playing impor- 
tant roles. For g = gc and T ^ the system approaches 
the quantum critical fixed point characterized by the ex- 
tinction of the spin- wave modes and the absence of long- 
range order. As shown in Fig. [21 the quantum critical 
region is defined by the crossover lines T ^ \g — gc\'^, 
where (j) = ^'^Si with v^ the 3d-Heisenberg correlation- 
length exponent. In particular, we find^" the following 
low-T behavior in the quantum critical region: 
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-l/z 
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f^ rpdj Z 
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with the standard result for the low-T quantum critical 
specific heat C and exponents satisfying scaling relations 
proper of this region^MS.. Similarly, for the T ^ 0, g < 




FIG. 2: Schematic RG {g,t oc T) flux diagram for the Id- 
quantum z = 2 NL a modeL Semiclassical {g* — T* = 0) and 
quantum-critical {g* = g^, T* = 0) fixed points are shown, 
as weU as the flux fines indicating their stabihty with respect 
to infrared perturbations. The segment T = Q,Q < go < 
gc, corresponds to the loci of points in which Lieb's theorem 
is included, with presence os stable ferrimagnetic states of 
quantum AB2 and alternate spin chains. 



gc semiclassical region, we find 
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where the specific heat is determined by the spin wave 
contribution. In Eqs. ([8]) and ([9]) the amplitudes of the 
observables cannot be completely fixed by the RG pro- 
cedure. We thus identify the asymptotic low-T critical 
behavior described by Eq. ^ with that of the quantum 
S — 1/2 AB2 and alternate spin ferrimagnetic chains, as 
well as that of the Id-quantum S = 1/2 Heisenberg fer- 
romagnet, such as the organic ferromagnetic compound 
p-NPNN*. In addition, an interesting question arises 
regarding the access of such AB2 chains in the half- 
filled strong-coupling limit to the quantum disordered 
and quantum critical regimes of the quantum z = 2 NL 
a model. This scenario, if accomplished, might involve 
the presence of extra frustrated couplings in the unit cell 
structure. In any case, we would like to mention that 
our predicted one-loop critical behaviors for the renor- 
malized classical and quantum critical fixed points are in 
agreement with those of the FM transition in Id itiner- 
ant electron systems in the context of a Luttinger liquid 
framework^^. However, while our localized spin disor- 
dered phase is gapped, the quantum disordered phase in 
Ref^ behaves as an ordinary gapless Luttinger liquid. 
Obviously, further theoretical work is needed in order 
to clarify the physical scenario predicted for the critical 
behavior of the quantum NL a model with a FM Wess- 
Zumino term due to the AB2 topology. 

In order to improve the understanding of the role of the 
quantum and thermal fluctuations, topology, and spin 
symmetry to the properties of the ferrimagnetic AB2 
chains, we have also performed a number of analytical 
studies using Ising, Heisenberg and spherical Hamiltoni- 
ans as model systems in this unit cell structure-^— ^. 



First, by regarding the spin operators in Eq. ([2]) as 
Ising variables, Sia = ifi/^, we apply^- the RG decima- 
tion of B sites and obtain the exact Gibbs free energy 
as function of the effective coupling J* and field H*. At 
zero field the ground state result J* = — 2J < for the 
effective coupling between A sites indicates the presence 
of a ferrimagnetic structure with average spins at B sites 
pointing opposite to those at A sites, implying in a unit 
cell average spin (S^^n) = h/2. As H increases at T = 0, 
we notice that J* increases linearly with H and vanishes 
for H > J; conversely, H* first decreases linearly with H 
for H < J, and then increases also linearly, changing sign 
at the critical field He = 2J. At H = He a. first order 
transition occurs with a discontinuous change of (S^^u) 
from h/2 to its saturated value 3h/2 (see Fig. ^. At fi- 
nite temperatures the described effects are less dramatic, 
and the unit cell average spin grows continuously with 
the field from at i? = (disordered state at finite T) to 
the saturated value 3h/2 as i? ^ cxd. The H — results 
are corroborated by the calculation of the two-spin cor- 
relation function, which is related to the susceptibility 
through the fluctuation-dissipation theorem. Actually, 
we have also found that X ^ ^ ^ sxT>[J/{kBT)]/{kBT), 
leading to the relation between the corresponding criti- 
cal exponents 'y = v = 2 — a, and from the behavior of 
the correlation function at T = and magnetization at 
T — and H -^ it follows that 77 = 1 and 6 — 00. This 
set of exponents belongs to the same class of universality 
of decorated Id ferromagnetic Ising systems^*. 

Now, by considering Heisenberg spins in a classical con- 
text, in which quantum fluctuations are absent, the low-T 
and high-T limits of the H = free energy and two-spin 
correlation functions have been calculated^. It is in- 
structive to compare the T — > result obtained for the 
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FIG. 3: Average spin per unit cell, {S^^n), in units of h, as 
function of the dimensionless magnetic field, H/ J , and di- 
mensionless temperature, t — ksT/J, for the ferrimagnetic 
S = 1/2 AB2 Ising chain. 



classical ferrimagnetic AB2 Heisenberg chain, 
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with that of the classical ferromagnetic linear Heisenberg 
chain^, 
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and that of Takahashi® for the quantum ferromagnetic 
spin-1/2 linear Heisenberg chain, 
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3C(l/2) {kBT^/^ , 3C'(l/2) keT 
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where C(l/2)/(27r)i/2 « -0.583. At low-T Fisher's and 
Takahashi's leading terms coincide and are three times 
larger than that of the AB2 chain, due to the unit-cell 
topology effect. Moreover, the second term in Eq. P^ . 
absent in the classical models, is related to the fixing 
of the anomalous entropy and specific heat classical be- 
haviors when quantum fluctuations are not present. We 
notice that this x ^ T~^ leading result as T ^ has 
also been obtained for the semiclassical fixed point of the 
quantum z = 2 NL a model, related to the quantum AB2 
Heisenberg chains (see above). 

By taking quantum fluctuations into account, we also 
calculate^ the three spin-wave modes of the quan- 
tum spin-1/2 AB2 Heisenberg model using the Holstein- 
Primmakov transformation and subsequent diagonaliza- 
tion via the Bogoliubov-Tyablikov method, namely, one 
non-dispersive optical, e^, one dispersive optical, e 
one acoustical mode, e^: 



^, and 



= J 
J 



+ H, (13) 

{±l + [l + 8sin^{ka)]^/^}TH ; 



notice in the acoustical mode the presence of a quadratic 
ferromagnetic dispersion relation e^ — 2J{ka)^, ka ^ 1, 
H = 0. From this result, corrections due to quantum 
fluctuations to the average values {S^) = —{S\) = 
h/2 are derived, although the result of Lieb's theo- 
rem, (Sceii) = fi^/2, remains true. In a mean-field ap- 
proach^i^, we relate the quantum thermal spin aver- 
ages at sites A and B to the respective Weiss molecular 
fields and find, in the simplest case in which they are 
assumed to be parallel to the z direction (Ising-like so- 
lution): (S"!;) = h/2 for all H > 0, and (S^) = -fi/2 
ioT < H < He, whereas {S%) = fi/2 for H > H^. 
We notice that He = 2J is the critical field below 
which the ferrimagnetic ordering is favoured, in agree- 
ment with the above result for the spin-1/2 AB2 Ising 
model. Indeed, the unit cell spin reads {S^^n) = fi/2 
for < i/ < He, and (S"^^,;) = 3fi/2 for H > He. On 
the other hand, in the case the x and y components are 
also considered, a quite interesting scenario emerges, with 




FIG. 4: Average spin at sites A (a), sites B (b) and per unit 
cell (c), in units of h, as function of the reduced field, H/Hc, 
for the quantum ferrimagnetic S = 1/2 AB2 Heisenberg chain, 
with He = 2J. (d) The field dependence of the Gibbs free 
energy shows that the continuous solution (solid line) for the 
magnetization is the stable phase. The Ising-like solution is 
shown for comparison (dashed line). 



(S-^) = :j^h/2 for < H < He/2, and (S^) = fi/2 for 
H > 3He/2, where the plus (minus) sign refers to a = _B 
(A) sites; for intermediate fields, He/2 < H < 3He/2, 
one has that {SX) = -h[3He/{4:H) - H/He]/2 and 
(5|) = h[3He/i8H) + H/{2He)]/2. These results im- 
ply in the unit cell average spin fi/2 for H < He/2, 
with ferrimagnetism sustained, and the saturated 3fi/2 
value for H > 3He/2, as in the Ising-like solution. 
A linear increase with H arises for intermediate fields: 
(S'eeii) = H/He, for He/2 < H < me/2 (see Fig. H. 
In this regime the average spin at sites A continuously 
rotates seeking a full alignment with H, accompanied by 
a rotation of the spins at sites B, such that the transver- 
sal spin components at sites A and B always cancel out. 
To achieve this cancellation the spins at sites B rotate 
in the opposite direction up to a maximum polar angle 
9 = 7r/6 and then rotate back [see Fig. Iljb)]. These re- 
sults are corroborated by the analysis of the Gibbs free 
energy [Fig. Sfd)] . 

At last, a spherical version of the quantum AB2 spin 
Hamiltonian has also been studied^. For this purpose, 
a chemical potential (/i) term is added to Eq. ^ in or- 
der to take care of the spherical constraint, ^iai'^fa) = 
N/A, where N is the total number of sites. Quantum 
fluctuations are introduced, associated with a quantum 
coupling parameter g, through a kinetic energy term 
(5/2) X^m -^iQ! i'^ which Pia are momentum operators 
canonically conjugated to each spin degree of freedom. 
By diagonalizing the Hamiltonian in a space of proper 
bosonic operators, we obtain two dispersive eigenmodes, 
also present in the linear AF spherical model, and a fiat- 



band induced by the AB2 topology. At the only critical 
point, g — T ~ H = Q, the ferriniagnetic long-range 
order is present with {Sb) — —{Sa)/V2 = fi.-\/3/4 and 
i'^ceii) — fii^~ l)-\/3/8. Interestingly, in the quantum 
AB2 spherical case the average spin per unit cell is less 
than H/2, in contrast with the result of Lieb's theorem 
for the AB2 Hubbard model in the strong-coupling half- 
filled limit and the quantum AB2 Heisenberg chain with 
AF couplings. Calculation of the correlation functions at 
g = T = H = show that they are distance independent 
and finite, consistently with the ferriniagnetic order and 
the spherical constraint. Outside this critical point, for 
any finite g, T or H, quantum and/or thermal fluctua- 
tions destroy the long-range order in the system, which 
in this case displays a finite maximum in the susceptibil- 
ity. In this regime spins remain ferrimagnetically short- 
range ordered to some extent in the {g, T, H} parameter 
space as a consequence of the AF interaction and the 
AB2 topology. Indeed, we notice in Fig.[5]that, although 
the field-induced unit cell average spin, {S^^n), displays 
quantum paramagnetic behavior for any finite g or T, the 
spins at sites A and B can display opposite orientations 
depending on the values of g, T and H. Therefore, for 
special regions of the parameter space {g, T, H} spins at 
sites A points antiparallel with respect to those at sites 
B, thus giving rise to a rapid increase in the unit cell 
average spin for very low H and a field-induced short- 
range ferrimagnetism, which is destroyed for large g, T 
or H. In addition, to better characterize the approach 
to the g = T — H — critical point, we have considered 
several paths. For T ^ and g — H = the suscepti- 
bility behaves as x ^ T^^, as also found in several clas- 
sical and quantum spherical and Heisenberg models (see 
above). On the other hand, for H ^ and g = T = 0, 
we find x ~ H~^, and for 5 — > and T — H — 0, 
X ^ exp(c(7^^'^), where c is a constant, evidencing an 
essential singularity due to quantum fluctuations. In any 
path, the relation x ~ ^^ is satisfied. We also mention 
that the known drawback of classical spherical models as 
T — > regarding the third law of thermodynamics (fi- 
nite specific heat and diverging entropy) is fixed in the 
presence of quantum fluctuations, g ^ 0. 



B. Numerical Results 

The ferriniagnetic ordering can be probed through the 
magnetic structure factor: 



'^(^) = ^E^'''"^'^'^(s..s,) 



(14) 



The condition for a long-range ferromagnetic ordering is 
that 5(0) ~ N; while ^(Tr) ~ A^ in a long-range an- 
tiferromagnetically ordered state; A ferriniagnetic long- 
range ordering fuUfill these two conditions: S{0) ^ N 
and 5(7r) ~ A^. This is the case for the AB2 Hubbard 
and Heisenberg chains as exemplified in Fig. [5] through 
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FIG. 5: (a) Average spin per unit cell, (Seen), and (b) spin 
averages at sites A (-0") 8.nd B (-□-), in units of h, as function 
of H/J, for g = 0.05J and T = 0.05J, calculated for the 
quantum spherical S = 1/2 AB2 model. Inset of (a): very- 
low-field regime. 



the exact diagonalization (ED) of finite clusters. Due to 
the critical nature of both chains at low temperatures, 
the correlation length ^ and x{q = 0) = S{q = 0)/{kBT) 
satisfy power law behavior: ^ ^ T~^ and x ~ T~'^ as 
T ^ 0. Since ^ ~ A^ at T = 0, using scaling arguments 
and the results of Fig. H we have T"'' - T-^/T, i. 
e., ^ — V = 1, in agreement with the values 7 = 2 and 
v = 1 derived using renormalization group techniques. 
Furthermore, the ferrimagnetism of the model was also 
manifested through Hartree-Fock and Quantum Monte 
Carlo methods'''^. 

Systems with a ferriniagnetic GS naturally have ferro- 
magnetic (lowering the GS spin) and antiferromagnetic 
(rising the GS spin) magnons as their elementary mag- 
netic excitations. The AB2 chain have three spin wave 
branches^: an antiferromagnetic mode (AF mode), de- 
fined as /^Es+{q) = Ei^S" = Sg + l,q) - Eas; and 
two ferromagnetic ones (Fl and F2 modes), derived from 
AEs-iq) = E{S' = Sg-\,q)- Eas, where Eas is 
the GS energy and E{S^,q) are lowest energies in the 
sector {S^,q}, with the lattice wave-vector q = 2Trl/Nc, 
where / = 0,1,..., Ac— 1. These modes are depicted in 
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FIG. 6: Magnetic structure factor S{q) for (a) the AB2 Hubbard chain with U = 2t and (b) the AB2 Heisenberg chain. The 
inset presents the size dependence of the ferromagnetic [5(0)] and antiferromagnetic [S(7r)] peaks. Dashed lines are guides for 
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Fig. [7] (a) for the Heisenberg model: the AF mode has 
a gap A5+ = 1.7591J; the gapless Fl mode is the Gold- 
stone mode, consistent with the symmetry broken phase 
of the chain; and the F2 mode has a gap As_. = 1.0004 J. 
The gapped F2 branch is flat and is associated with the 
formation of a singlet state between the B sites in one 
cell, while the other cells have B sites in triplet states, 
as illustrated in Fig. [7] (b). The localized nature of the 
excitation is associated with the Hamiltonian invariance 
under the exchange of the B sites of any cell, this sym- 
metry implies that the many-body wave function has a 
definite parity under the exchange of the spacial variables 
associated with these sites. Since these dispersive modes 
preserve the local triplet bond, they are identical to those 
found in the spin- ^ /spin- 1 chain™. Surprisingly, Linear 
Spin Wave Theory^ predicts that A5_ — 1, very close to 
our estimated value: As- — 1.0004J. Moreover, a good 
agreement is found for the gapless Fl branch in the long 
wave-length limit. However, both LSWT and mean field 
theory^ predicts A5_|_ = 1, deviating from our estimated 
exact diagonalization value: As_|_ = 1.7591J, which is in 
excellent agreement with numerical and analytical calcu- 
lations for the spin- i /spin- 1 chainl^. On the other hand, 
the Interacting Spin Wave Theory^ derives a better re- 
sult for As_|_, but it implies in a higher shift for Ag^ (flat 
mode) not observed in our data of Fig. [7| (a). 

On the other hand, the AF mode is relevant in the anal- 
ysis of the response to an applied magnetic field H. The 
AF gap found above is responsible for a plateau in the 
curve of the magnetization per spin [m{H) = {S^)/{Nh)] 
as a function of H. In fact, it has been shown- that if 
j/(s — m) = integer, a plateau may appear in the mag- 
netization curve of the Heisenberg model. In the last 
equation, s is the site spin quantum number and v is the 
number of sites in one unit cell of the GS for a given value 
of H. The AB2 Heisenberg chain has s = 1/2 and three 
sites per unit cell {v = 3); so, unless the system spon- 
taneously breaks the translation invariance, we expect 
plateaus at m = 1/6 and m — 1/2. This is indeed what 
is observed in Fig. [H) The plateau width at ?ti = 1/6 is 



exactly given by A5+ , and is a measure of the stability of 
the ferrimagnetic phase. For higher fields, the magnetiza- 
tion increases in the expected way-^, as shown by the full 
line in Fig. [H before saturation at m = 1/2 for H = 3J. 
This field-dependent behavior contrasts with the linear 
one predicted by mean-field theory^^ shown in Fig. [31 
Exact diagonalization results''^ indicate that the antifer- 
romagnetic spin gap and, consequently, the plateau at 
m = 1/6 exists for any finite value of U, with the plateau 
width (As+) nullifying as C/^ in the limit C/ — > 0. 

Away from half-fiUing^i^, the AB2 Hubbard model 
exhibits a rich phase diagram depending on the electronic 
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FIG. 7: (a) Ferromagnetic (F) and antiferromagnetic (AF) 
spin wave modes of the Heisenberg AB2 chain for Nc = 10 
(circles), 8 (triangles down), 6 (triangles up). Solid lines are 
the Linear Spin Wave results from Ref.— ; dashed lines are 
guides to the eye. (b) Illustration of the F2 mode: ellipse 
indicates a localized singlet state. 
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FIG. 8: Magnetization as a function of the dimensionless ap- 
plied magnetic field H/J for the AB2 Heisenberg chain. The 
full line is a curve traced from the midpoints of the steps 
found in the finite size results, except at plateau regions. 



density, n = Ne/N , or doping 5{— 1 — n) from half-filling, 
and the Coulomb coupling U . The doped region was ana- 
lyzed through Hartree-Fock^-, exact diagonalization^'^-^^ 
and density matrix renormalization group (DMRG)^^, 
which is the state-of-the-art method^^ for the study of 
the GS properties of one-dimensional quantum lattice 
models. DMRG results suggest that in the undcrdoped 
region and for U — 2i, the ferrimagnetic phase sustains 
up to ^ ~ 0.02 while for 0.02 < 5 < 0.07 hole itiner- 
ancy promotes incommensurate spin correlations (a spi- 
ral phase^^) with a 5-dependent peak position in the mag- 
netic structure factor, as shown in Fig. [§] (a) and (b). 
For U ~ 00 and 5 — Q the GS total spin is degener- 
ate, whereas for < i5 < 0.225 hole itinerancy (Nagaoka 
mechanism2iiZ^) sets a fully polarized GS, as shown in 
Fig. [SI (c). 

For higher doping, the system phase separates^ into 
coexisting metallic and insulating phases for Sps(U) < 
S < 1/3 [with 6ps{oo) w 0.225 and Sps{2) w 0.07]. In 
fact, the Hartree-Fock solution is unstable in this region 
and a Maxwell construction is needed^'^. The local par- 
ity symmetry is even (odd) in the insulating (metallic) 
phase. In Fig. [TU] we present spin correlation functions 
at S = 0.18 {6 = 0.28) and U = 2 {U = 00) calculated 
through DMRG (for these parameter values, the system 
is found in a phase separated state). In Fig. [TOKa) the 
correlation function between the spin at the extrem site of 
the phase with odd parities and the others spins, (Si -S;), 
evidences the spiral phase for U = 2t. For U = 00 the 
local magnetization at sites A and Bi + B2 shown in 
Fig. [TUTb) clearly displays the coexistence between the 
Nagaoka ferromagnetic phase and a paramagnetic one. 
Notice that, in the Nagaoka phase, the magnetization dis- 
plays spatially modulated profiles due to hole itinerancy. 
Further, the local parity symmetries of the two coexisting 
phases are manifested in the correlation function between 
B spins at the same cell, (S^^ -Ss^), shown in Figs. \Wi c) 
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FIG. 9: (a) and (b): Magnetic structure factor for U — 2t 
and A'" = 100, using DMRG, in the underdoped region, (c) 
Total spin per cell Sg/Nc as function of S for 17 = 00; the 
variational approach is described in Ref.— 



and (d) : in the phase with odd (even) parities cell triplet 
(singlet) states predominate. In Fig. [TUf e) we illustrate 
the phase separation for U = 00. DMRG and exact di- 
agonalization results^S. indicate that the phase separation 
region ends precisely at 5 — 1/3. For this doping the elec- 
tronic system presents finite spin and charge gaps with 
very short ranged correlations and is well described by 
a short-ranged resonating-valence bond (RVB) state—, 
with the electrons correlated basically within a cell, as il- 
lustrated in Fig. [TUTf). A crossover region is observed for 
1/3 < 6 < 2/3, while a Luttinger-liquid behavior^* can 
be explicitly characterized for S > 2/3. Luttinger liq- 
uids are paramagnetic metals in one-dimension exhibit- 
ing power-law decay of the charge and spin correlation 
functions and the separation of the charge and spin exci- 
tation modes. In particular, the asymptotic behavior of 
the spin correlation function is given by 



CllH) 



cos{2kFl)[Hl)]^/'^ 



(15) 



where kp is the Fermi wave vector and Kp is the model- 



dependent exponent. The predicted behavior in Eq. (fT5|) 
fits very weU the DMRG data at S = 88/106 using 
Kp = 0.89 {U = 2t) and Kp = 0.57 [U = 00), as shown 
in Fig. [11] We remark that these values are close to 1 
(noninteracting fermions) and 1/2 (noninteracting spin- 
less fermions) for U = 2t and U = 00, respectively. 

In addition, we mention that the commensurate doping 
S — 2/3 is insulating, with a charge gap nullifying with U 
in a similar manner as the one of the Hubbard model in 
a linear chain at half-filling2^, while the spin excitation 
is gapless. 
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IV. CONCLUSIONS 



In this contribution on the celebration of the 80th 
birthday anniversary of Prof. Ricardo Ferreira, we have 
presented a brief review of the main experimental and 
theoretical achievements on quasi-one-dimensional mag- 
netic compounds, featuring those with AB2 unit cell 
structure. 

As reported, this has been an area of intense activity, 
particularly since the first experimental announcements 
in the mid 80 's'^^'^'^. Nowadays several groups all over 
the world, involving chemists, physicists, and material 
scientists, are engaged in the characterization and de- 
scription of properties of the already known materials, 
as well as doing great efforts towards the design and syn- 
thesis of new compounds, with novel properties suitable 
for technological applications^. 

From the fundamental point of view, these compounds 
have been used as a laboratory in which many theoretical 
concepts and predictions in the field of low-dimensional 
materials have been tested. 

In conclusion, it seems clear that this interdisciplinary 
research area will remain an exciting and topical one for 
many years to come, offering new challenges both from 
the scientific and technological aspects. 
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